#轮函数F
def F(x,y):
    u=xor32(x,y)
    u1=u[0:8]
    u2=u[8:16]
    u3=u[16:24]
    u4=u[24:32]
    S={'0':5,'1':10,'2':15,'3':4,'4':9,'5':14,'6':11,'7':8,'8':2,'9':7,'10':12,'11':13,'12':3,'13':6,'14':1,'15':0}
    v1=format(S[str(int(u1[0:4],2))],'04b')+format(S[str(int(u1[4:8],2))],'04b')
    v2=format(S[str(int(u2[0:4],2))],'04b')+format(S[str(int(u2[4:8],2))],'04b')
    v3=format(S[str(int(u3[0:4],2))],'04b')+format(S[str(int(u3[4:8],2))],'04b')
    v4=format(S[str(int(u4[0:4],2))],'04b')+format(S[str(int(u4[4:8],2))],'04b')
    v=v1+v2+v3+v4
    z=xor32(xor32(v,demo(v,3)),demo(v,10))
    return z
